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We introduce matching functions as a means of summing heavy-quark logarithms to any order. 

Our analysis is based on Witten’s approach, where heavy quarks are decoupled one at a time 
in a mass-independent renormalization scheme. The outcome is a generalization of the matching 
conditions of Bernreuther and Wetzel: we show how to derive closed formulas for summed logarithms 
to any order, and present explicit expressions for leading order (LO) and next-to-leading order 
(NLO) contributions. The decoupling of heavy particles in theories lacking asymptotic freedom is 
also considered. 
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I. INTRODUCTION 

Decoupling a heavy quark when the renormalization 
scheme is mass independent was originally discussed by 
Witten |l||. He showed that the results can be elegantly 
expressed in terms of a renormalization group (RG) in¬ 
variant running coupling cth associated with the mass 
rrih of the heavy quark h. Subsequently, Bernreuther and 
Wetzel [2, Li [1IH proposed a systematic method for deal¬ 
ing with the matching problem, i.e. the lack of explicit 
decoupling in mass independent schemes. They applied 
the Appelquist-Carrazone decoupling theorem |(| to the 
gluon coupling Oq° in the momentum subtraction (MO) 
scheme, i.e. renormalized at space-like momentum Q 

a Q°Lith h = a Q°L h + °^ m h 1 ) (!) 

and compared calculations of ag° in the full F = f + 1 
and effective / flavor MS (modified minimal subtraction) 
theories. When C^m^ 1 ) terms are neglected, the strong 
coupling otp = <?|,/(47r) for the E-flavor MS theory is 
calculable as a power series in its /-flavor counterpart af 
and logarithms of rrih- Results for the first few loops of 
perturbation theory appear in the literature nidi 
0, la, 0 • Bernreuther has constructed a similar matching 
procedure to deal with the effects of mass renormalization 
1 

This paper arises from the observation HQ that the 
RG relates coefficients of perturbative mass logarithms 
~ a r F ln s rrih in matching relations. This suggests that 
we seek an analogy with the behavior of Green’s func¬ 
tions at large momenta q , where in general EMI , each 
perturbative order in the Gell-Mann-Low function T (x) 
or the Callan-Symanzik /?, 7, S functions determines all 
coefficients to logarithmic order k, i.e. all coefficients of 
momentum logarithms {a(p +K ln r q, r = 0,1, 2 ,...}. 

We set up the formalism for mass logarithms by intro¬ 
ducing matching functions F(a.h) and Q(cth) associated 


with coupling constant and mass renormalization. The 
coupling constant and masses are matched to all orders, 
with heavy quarks (in this article) decoupled one by one 
(.F = f + 1). For coupling constant renormalization, we 
derive in Sec. rm the following key equations, 




1 

(3f(x) 




1 

/3 f (x) 


+ FF^f(ah) 


( 2 ) 

( 3 ) 


where rrih is Witten’s renormalization group invariant 
heavy-quark mass pj], (3f(x) is the MS /3-function for the 
/-flavor theory, and 

fi = MdimV / 47re -7 ^ 2 , 7 = 0.5772 ... (4) 


is the MS scale derived from the scale puim used to de¬ 
fine dimensional regularization and renormalization. The 
matching function J-F^f is a series in oih whose coeffi¬ 
cients can be determined perturbatively by comparison 
with ©. The desired matching relation between ap and 
af is the result of eliminating ah from Eqs. 0 and 0: 

a F = a F (af,\n(rh h /fl)). (5) 

Similar conditions for mass matching are presented in 

Sec. ivm 

The role played by T, Q and the /3,7, 5 functions in 
matching conditions is just like that of the /3,7, 6 func¬ 
tions for large-momentum logarithms. Each order of per¬ 
turbation theory for these functions determines the co¬ 
efficients of a new tower of mass logarithms: leading or¬ 
der (LO), next-to-leading order (NLO), next-to-next-to- 
leading order (NNLO), and so on. Both F(dth) and Q(ah) 
vanish for LO and NLO, but then there are contributions 
from successive terms in their power series in ah, starting 
with NNLO for F(oth) and NNNLO for CJ(ah)- 
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We find, as for large-g logarithms, that results for cou¬ 
pling constants, running couplings and light masses can 
be most elegantly expressed as closed expressions or gen¬ 
erating functions for summed towers of mass logarithms. 
Examples are the decoupling formulas for a/+i and ah 
quoted previously by us [l2|, which we derive in Sec. eh 
Almost all of our results are for quantum chromodynam¬ 
ics (QCD) with three colors, but the technique can be 
applied to any renormalizable theory. 

Sections UH and rrm are brief summaries of Witten’s 
treatment of heavy-quark decoupling in QCD and the 
matching procedure of Bernreuther and Wetzel for cou¬ 
pling constants. This lays the foundation for the RG 
analysis in Sec. E3 from which we are led to construct the 
matching function T for coupling constant renormaliza¬ 
tion. Perturbation theory for T is considered in Sec. El 
with the result that the first non-zero term (NNLO) in 
T is obtained. In Sec. ED we show that Eqs. © and 
© lead directly to closed expressions for heavy-quark 
logarithms to a given logarithmic order, and present 
explicit NLO expressions. Section EDI is an extension 
of our RG analysis to deal with the matching problem 
for mass renormalization. It is here that we introduce 
the mass-matching function Q. In Sec. IVIIII we decou¬ 
ple more than one heavy quark sequentially, for exam¬ 
ple In (rrit/p) In (nib/ft) —> oo, and derive the NLO 
closed formula for coupling constant renormalization in 
this limit. 

Section [K] suggests that the consistency of theories 
lacking asymptotic freedom, such as quantum electrody¬ 
namics (QED) with heavy leptons, be tested by imposing 
the physical requirement that all heavy particles decouple 
in the infinite-mass limit. Both ultra-violet and infra-red 
stable fixed points enter the analysis. 

Other applications of our technique are discussed in 
the concluding Sec. El 

II. WITTEN’S METHOD 

This section summarizes some key points of Witten’s 
procedure 0112 

By convention, the same MS scale /i is used for the ini¬ 
tial E-flavor and all residual /-flavor theories. Whenever 
heavy quarks (masses to/j) are decoupled, 

E —> / flavors, tof —> oo 

all parameters of the residual /-flavor theory are held 
fixed: the scale / 2 , all momenta p, the coupling a/, and 
all light-quark masses mif. In any order of perturbation 
theory, amplitudes 

Af = -4f(p, P, a.Fi miF, m h ) 

are power series in to / 1 with each power modified by a 
polynomial in \n[rrih / p) . We will consider the leading 
power Af- 

Af = Af{ 1 + O^m'p 1 )}. (6) 


The notation 0{rn h l ) refers to any sub-leading power, 
including its logarithmic modifications. 

Logarithms in Af for rrih ~ oo are generated by 1PI 
(one-particle irreducible) subgraphs with at least one 
heavy-quark propagator and with degree of divergence 
at least logarithmic. It is as if all contributing 1PI parts 
were shrunk to a point. All F- flavor amplitudes Af tend 
to amplitudes Af of the residual /-flavor theory, apart 
from TO^-dependent renormalizations of the coupling con¬ 
stant, light masses, and amplitudes [f|: 

A.F(p,ft,aF,miF,m h ) 

= ^2z AA '(a Fl m h /p)A l f (p,p,af,mif) ( 7 ) 

A' 

af = af(a F ,m h /p) , mif = mi F D{a F ,mh/p )• (8) 

For practical applications, Eq. © has to be inverted, so 
that af and to;/ become the dependent variables instead 
of Of and to/f- That is because we hold a/ and to;/ 
fixed as mh —> oo. 

For any number of flavors / (including E), let 

Vl = ii k + Ma,) ^ +S ' {a (9) 

be the corresponding Callan-Symanzik operator. Since 
Af satisfies an E-flavor improved Callan-Symanzik equa¬ 
tion 0 , so also does its leading power: 

{F>f + 7f(o:f)}-^f = 0- (10) 

In general, both 7 p and Z = (Z AA >'j are matrices. 

If we substitute © in m and change variables, 

V F = p^+{V F a f )^ + ^{V F m k f )^- ( 11 ) 

the result is an improved Callan-Symanzik equation for 


each residual amplitude, 

{V f +^f(af)}Af = 0 (12) 

where the functions 0 i 

Pf( a f) = F> F a f (13) 

^/( a /) = Fplnmif (14) 

7/(a/) = Z- 1 ( lF {a F ) + V F )Z (15) 


depend solely on af. The absence of to; dependence in 
the renormalization factors in Q and © ensures mass- 
independent renormalization for the residual theory. 

While these equations hold for any / < F, their so¬ 
lutions can be readily formulated in terms of running 
couplings only when the heavy quarks are decoupled one 
at a time. Indeed, Witten’s running coupling 

a h = a h (a F ,ln(m h /p)) 


(16) 
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is defined for the case F = / +1 where just one quark h is 
heavy, with MSi? renormalized mass mh- The definition 
of ah is formulated implicitly Jj]: 

r&h. 

In (m h /p,)= / dx (l - S F (x)) / fi F {x). (17) 

J Oip 

It satisfies the constraints 


where 


r 


h —loop 


— [ dss(l — s) In 

Jo 


s(l-s)Q 2 

ft 2 


( 22 ) 


is the contribution of the heavy-quark loop, and r ot her 
represents other one-loop terms. 

The leading power contributed by the heavy-quark 
loop is 


a h (a F ,0)=a F , a h (a F ,oo)=0 (18) 

where the latter follows from the asymptotic freedom of 
the F -flavor theory (F < 16). Eqs. 0, (THil) and m 
imply that ah is renormalization group (RG) invariant: 


T/i-ioop = Clo In —r- 4- Cnlo + 0(Q 2 /to 2 ) (23) 

with coefficients for leading and non-leading logarithmic 
orders given by 


V F ah = 0. 


(19) 


Clo = 1/(37t) and Cnlo = 0 . (24) 


III. MATCHING COUPLING CONSTANTS 

Generally, the solutions of Witten’s equations depend 
on renormalized parameters a F and mi F of the original 
F = f + 1 flavor theory, whereas the limit mh —> oo is to 
be taken with parameters a / and mi / of the residual the¬ 
ory held fixed. To complete the analysis, it is necessary 
to derive asymptotic series in ln(?nft//2) which relate the 
initial and residual parameters, i.e. to “match” a F and 
mi F with af and rnif- As noted in Sec. 3 Bernreuther 
and Wetzel 3,(3, 401 have set up a systematic procedure 
for this. This section is a brief account of their scheme 
for the case of coupling-constant matching. 

The decoupling formula 0 works to any order of per¬ 
turbation theory, so the task is to express the leading 
power of the RG-invariant gluon coupling ag° with and 
without the heavy-quark h as perturbative series in a? 
and af respectively. Generally this involves gluon and 
other self-energy insertions and a vertex amplitude such 
as fermion-gluon @0 or ghost-gluon d|. 

For one-loop contributions [14(, vertex and propaga¬ 
tor corrections cancel {Z\ = Z 2 ), so only the gluon self¬ 
energy amplitude 


The vanishing of the NLO constant term is a well-known 
characteristic of the MS gluon self-energy Ill- 

Eliminating r o ther from Eq. <EU and combining the 
result with Eqs. {HI and {23, we recover the standard 
one-loop matching condition 

a f+i ~ a f 1 = ln + °( a fi 1 ) (25) 

or equivalently 

2 

a f+i = a f - + °( a % m h 1 )- ( 26 ) 


The two-loop analysis is much more complicated, so 
we simply quote the result 0 S 0,0, taking account of 
a subsequent correction 0, 0 fej. We find it convenient 
to consider the inverse form where cn/+i is written as 
a series in af. For the special case of three colors, the 
result is: 


a f 1 m h a f 


hr 


lla 3 , 


/ 


r>' l 11 

24-7T 2 fi 


77 if lla 3 , . 


is needed. In that case, we can make the replacement 

“0°l, n™ —* “//(! - n /(<3» 

in Eq. 0 , with the result 

a f +1 — ot f+1 H f+1 (Q) 

= af 1 -af 1 H f (Q)+0(mh\a 2 f ). ( 20 ) 

Comparing the original and residual theories, we have 


The first three terms of the right-hand side belong to 
the leading order LO, i.e. they are proportional to af 
times a power of {af In(m/,//i)}. Only the fourth term is 
NLO; there is no 0(a 2 ) term independent of mh because 
the NLO constant in Eq. {23 vanishes. The fifth term 
is 0(a 3 f) and TOft,-independent, so it is the first example 
of a NNLO term. The three-loop result, including the 
NNNLO constant term, is now known |gj. 

Now we would like to know what the renormalization 
group implies for matching relations of this type. Some 
results for coefficients to a given order of perturbation 
theory already appear in J2, 3, 0 ■ Consider Eq. (3) of 
Ref. 0 , 


II/+l(f)?) — {T/j—ioop A Tother} 4“ G(oj_[_i) 

Hf{Q)=a f r othei + 0(a 2 f ) (21) 


aj_ _ a f+ 1 
7r 7r 


Wa/+i\ 


k =1 


Cfc (l n ^ L ) +0 ( m ^ 1 ) ( 28 ) 








4 


where Ck is a polynomial of degree k , as noted below 
Eq. (8) of Ref. 0 :' 

Ck = Ck,k (ln(m^//x 2 )) + ••• + Cfc, o- (29) 

The constants cqo, C2,o, 03,0, ■ ■ • are the remainders left 
when all terms depending on In(m 2 //x 2 ) are subtracted 
from the leading-power functions C\, C2, C3,.... Then, 
if all coefficients and RG functions are known to k — 1 
loops, the RG determines all /c-loop coefficients Ckj in 
Ck except for Ck,o■ The latter cannot be deduced from 
the RG, to any number of loops; rather, Ck,o must be 
calculated explicitly via a separate fc-loop matching cal¬ 
culation. For example, the NNLO coefficient — ll/( 727 r 2 ) 
in m is just —C2,0- 

Instead of Eq. ( 1251 ) . we prefer to consider the inverse 
relation 

OO 

a f+i = +J2 ak f +lpk ( ln ~= 1 ) +°( m h 1 ) ( 3 °) 

k= 1 ^ 

because that is what is required in order to take m-h —> 00 
with af held fixed. The analogue of Eq. ( 1291 ) is 

Pk = Pk,k (In (mh/fi)) k + Pk,k -1 (in {m h /p)) k 1 

+ • • • + Pk, 0- (31) 

An analysis in the style of Bernreuther and Wetzel pro¬ 
duces identical conclusions for the remainder constants 
p k ,o- given pi,o,P2,o, • ■ -Pk-i,o, one can use the RG to 
deduce p k ,k , ■ ■ -,Pk, 1 but not p k ,o- 

Most practical applications require that terms of the 
same logarithmic order be summed. This is straightfor¬ 
ward for LO logarithms, because the LO coefficients Ck,k 
in lt 251 ) obey a simple relationship 0 


Beyond LO, formulas for all the relevant coefficients 
become complicated, making order-by-order summation 
too cumbersome to be practical. The rest of this paper is 
concerned with a RG analysis which allows us to consider 
matching relations to a given logarithmic order without 
having to expand in perturbative order. 


IV. MATCHING FUNCTION 

Any RG analysis of decoupling involves ar least two 
renormalization groups: one for the initial E-flavor the¬ 
ory, and one for each /-flavor theory produced as a heavy 
particle decouples. We append a flavor subscript to make 
the distinction, viz. RGf or RG/. 

A key observation is that any quantity which is RGf 
invariant must also be RG/ invariant (/ < E). For exam¬ 
ple, Witten’s RGf invariant running coupling 5/ must 
satisfy the condition 

V,S h = ^+/3^S h =0. (37) 

Generally, the substitution 

V F —* V f (38) 

works when applied to any quantity which survives the 
limit irih —> 00 . An example is the formula 

2^60.5 = 2 / 50:5 = , 85 ( 05 ). (39) 


c k ,k = (ci,i) fe (32) 

which makes the series geometric: 

^£ = o“' +i /( 1 '^r ln T)' (33) 

This expression is leading order (LO) with respect to 
powers of o/+i = a F and \n(mh/p )■ Eq. (1331) implies 
that the term 0(o/, mj) 1 ) in (12511 is NLO or higher order: 


-1 


-1 


m h 


o r 1 1 o r — In 
' ' LO 3?r ll 


(34) 


This leads directly to the inverse of (I551i . 


«/+> ^“'/( 1 + i ( 1 " if ) 


LO 


(35) 


which agrees with the general result G3- However, the 
converse is not generally true. For example, the top- 
quark mass nit is RG / = 5 invariant, but it is certainly not 
RG / = 6 invariant. Therefore a study of the RG for the 
original E-flavour theory is both necessary and sufficient 
for the full implications of the RG to be understood. 

Our starting point is Witten’s definition of the 
invariant running coupling 5^. Let us regard this as a 
formula for In (mh/fi) in terms of and a F . Specifi¬ 
cally, the right-hand side is an integral from a F to clh 
involving RGf functions (3 F and S F . Can a similar for¬ 
mula be constructed from RG/ functions such that this 
mass logarithm becomes a function of ah and a/? 

If such a formula exists, it must be consistent with the 
requirements of the RGf group for the original theory. 
However, mass renormalization produces an unwelcome 
E dependence in equations such as 


where now LO refers to powers of a/ and In (mh/p). Note 

that the LO coefficients pk,k in Eq. (13111 are given by Ef ln(mft,/ p.) = 6 F (a F ) — 1. (40) 


Pk,k = (—1 / 37r) ^. 


(36) So let us amend the proposal: instead of the MSf mass 
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TOft, consider Witten’s invariant mass 1 


rah 

m h =m h ex p/ dx S F (x)//3 F (x). (41) 

J OLF 

Since fhh is RGf invariant, 


Vpirih = 0 (42) 

replacing nih by fhh in (TO eliminates the unwanted de¬ 
pendence on Sp: 


V F ln(m/i//i) = -1. (43) 

Notice that the formula m for \n(rrih/fi) is an immediate 
consequence of the definitions CD and CD of ah and 
fhh- As a check, Vp can be applied to the right-hand side 
of |2l to give the result —1, in agreement with Eq. (E31l . 

Now observe that, because of Eq. d. the replacement 
F —> / everywhere on the right-hand side of m produces 
a quantity which transforms in the same way under the 
RG f of the original theory: 


V F 



1 


= - 1 . 


(44) 


Comparing Eqs. m and d , we see that a RG f invari¬ 
ant quantity T can be defined as follows: 

fhh f^ h 1 

In - = - = / dx —— + T => V f T = 0. (45) 

P Ja f Pf{x) 

Since T is dimensionless, it can depend on ah-, but 
RGf invariance forbids dependence on other dimension¬ 
less variables, such as a F , a/, fhh/mh , or m-h/p- We call 
it the matching function: 


T = Fp^fiah). (46) 

This completes the derivation of the key equations m 
and 01: 


. rn h 

in 

fhh 


(•a-h 

dx 


Pf(x) 


In 


fc 

P Ja s 


a h 

dx 


P f {x) 


+ Fp^f(ah)- 


When ah is eliminated from these equations, the desired 
matching relation 0 is obtained: 

a F = a F (af,\n(fh h /p)), F = f + 1. 


1 See Eq. (16) of Q. Similar effective masses have been invented 
for the cases of large momenta El and light quarks Il5l. Their 
RG invariance mak es t hem useful in phenomenology pi 1 15(1 and 
lattice calculations G2. 


V. PERTURBATION THEORY FOR THE 
MATCHING FUNCTION 


Perturbative matching relations can be used to deter¬ 
mine successive coefficients in the Taylor series of the 
matching function FF^f(ah)- One needs to calculate 
the mass-independent terms Ck,o in lt29l) or equivalently 
Pk, o in <ed. As noted in Sec. rrm these constants cannot 
be deduced from the RG. This corresponds to the fact 
that cannot be deduced from the RG because it is 

RGf invariant. 

To illustrate the procedure, let us deduce the conse¬ 
quences for of the perturbative matching rela¬ 

tion ED- We need the two-loop /3-function for three 
colors: 

Pf(x) = -^(33-2/)- I |^(153-19/) + 0(x 4 ). (47) 
Its reciprocal is 


{Pfi x )} 1 


6n /I 
33-2/ 



where bf stands for the constant 


, 1 153-19/ 

bf =-— 

1 2tt 33 - 2/ 


O(x) (48) 


(49) 


and b'j is another constant whose precise value is not of 
concern here. 

The expansion (TTSll inserted into 0 and 0 yields the 
following equations: 


33 — 2(/ + 1) 

67T 


33-2/ 

67T 


. nih —i —1.7 , a h 

\n— = a h -a f+1 + b f+1 ln- - 

p J &f+1 

+ 0(a 2 f ) (50) 

In = aT 1 — af 1 + br In — + 0(a 2 f ) 
H 1 a f 

+ —-p——Ff + I^ s (a h ). (51) 


Note that the corrections are 0(a 2 ): contributions to 
m and CD from the constant term in (TIHl) are 

b'f +1 (ah - a/+ i) = 0(a 2 f ) and b' f (a h - a/) = 0(a 2 ) 


because ah-, a /+1 and af differ by 0(a 2 ). 

The next step is to eliminate cth- This is partially 
achieved by subtracting m from (ED: 


1 , rnh _i -i , . oih . , ah 

— in — = a u ,- a, + Ofin -6f+iin- 

37t p IX 1 af a f+i 

-\ - ——Ff +1 ^f{a h ) + 0(a 2 ). (52) 


Since \n{ah/af ) and \n(ah/otf+i) are both 0(a/), Eqs. 
ED and ED imply 


Ff+i^f( a h) = 0(a f ) 
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and so, from Eqs. (0 and 0 , we conclude 

a /+i / a h = 1 + ~ ^33 — 2 (/ + l)^j In ——b 0(a 2 ) 

a f /a h = l + ^(33-2/)ln^ + 0(a^). (53) 


The logarithms of these expressions can then be substi¬ 
tuted back into Eq. 0 , with the result 

1 , fn h _i _i 19a/ m h 
— ln — = “/+!-“/ — 


37T fi 




33-2/ 

67T 


JP /+w (5 /i ) + 0(a)). (54) 


The next step is to relate the logarithms of mh and 
mh■ First substitute the three-color formula 


S f( x ) = -—+0{x 2 ) 

7T 

into the definition m of mh, 

12 , a h 


, mh 

in- = 


In ——-b 0{a.f) 


(55) 


(56) 


m h 33 - 2(/ + 1) a/ + i 

and then substitute Eq. (1531 for cif+\/oLh- The result is: 


ln - 


m h f, 2a f \ 


7r / 


\n^ + 0(a}). 


(57) 


Then the logarithm of mh can be eliminated from 
Eqs. 0 and 0 : 




a' 


a ' +1 =a '“S 1 "T + S? 1 " 


+ j/(33-2/)J-, + 1 _ / (Si)+0(o}). (58) 

Comparing this with the two-loop matching condition 
0 , we see that all mass logarithms are correctly repro¬ 
duced, and that the first non-zero term in the matching 
function can be deduced from the constant NNLO term 
in (12711 : 


/i 




llQ / m h 

12tt 2 a 


Ff + i^f{a h ) = - 


11 


12tt(33 - 2/) 


a h + 0(a 2 ). (59) 


The 0(af l ) term in T can be found by substituting 
(0 back into Eq. 0) and repeating the above process 
using the three-loop /?/ and two-loop Sf functions. The 
answer follows by comparison with the known three-loop 
matching condition Q. That is the limit of current calcu¬ 
lations, but in principle, this strategy could be pursued to 
any order, with all mass logarithms correctly reproduced. 


VI. CLOSED EXPRESSIONS FOR 
HEAVY-QUARK LOGARITHMS 


The importance of the matching function T is that it 
allows us to work to a given logarithmic order without 


having to sum mass logarithms order-by-order in pertur¬ 
bation theory. Indeed, the role of T is essentially the 
same as that of the RG functions f3, 7 and 5: each term 
in the series for T corresponds to a particular logarith¬ 
mic order. For LO and NLO, T does not contribute, but 
NNLO requires that the O(a^) term in (I5H1) be included, 
NNNLO requires the O(a^) term, and so on. 

To illustrate, let us derive the closed NLO formula for 
the matching relation between a /+1 and a/ which we 
announced in [l 2 | . 

As in the previous section, we insert the expansion 0 
into © and m, but this time we omit the NNLO term 
Ff+i^r- 


33 — 2 (/ + 1) fhh 

p, nlo 


67T 


In— - a h 1 -a f l 1 + b f+ 1 


ln 


Oth 


a f +1 
(60) 

33 — 2/ fhh —— x —1.71 a h 

■in- = a h —a f + b f ln —. (61 

jjL NLO n 1 Oif 


67T 


The difference of these two equations is 

1 , mh _i _i ,, 7 , i a/ 

— in- = 07 ,,—a, — (Of — bf+i) In — 

3 ?t /i nlo W f w J+ ’ a h 


- bf+i ln 


a f 

a f +1 


(62) 


The logarithms on the right-hand side are already NLO, 
so we can use the LO parts of 0 and 0 to approxi¬ 
mate their arguments: 

otf/a h = 1 + ^(33 - 2/)ln^ (63) 

LO 67T /t 

a f/ a f +i + ( 64 ) 

LO o7T fl 

The result is a NLO generalisation of (1351) : 


a /+i 


NLO 


/ r af m h , , , , a/ m h 

“'/i 1 + sf ln + “'*'+>H 1 + u ln x 


'i 


37T /i 
+ a f(bf-b f+ i)!n 


37t /i 7 

l + ^(33-2/)ln^ 

07r n 


]}■ 


(65) 


If desired, \n(fhh/p) can be eliminated in favor of 
ln (mh/jj-)- The leading NLO effects of mass renormaliza¬ 
tion are due to the one-loop term of <5/ given by Eq. (17^1) . 
When this term is substituted into the definition (1411) of 
fhh , keeping all logarithms to this order, we find an ex¬ 
pression 


, m h . m h 

in- = in- 


- _ + Tj- 57 ( ln - ln ( 66 ) 

/i nlo /i 31 — 2/ V a/+i ah > 

to which Eqs. 0 and 0 can be readily applied: 


ln 


m h 


= ln^ + 


12 


fj, NLO /i 


, , °if , m h \ 

■H 1 + 


12 


31-2/ 


31-2/ 

ln(l + g(33-2/)ln^). (67) 
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So, by combining ED and ED, we arrive at a complete 
NLO formula for the matching condition: 


OL /-ii — OL f 
/ + i NLO J 


1 , m h 
— hr — 
3tt p 


c f 


In 


, i m h 

1 + -Mu — 
3tt p 


c f = 


+ df In 

142 - 19/ 
2tt(31 —2/) ’ 


l + ^(33-2/)ln^ 

07T fi J 


d f = 


57+16/ 


2tt(33-2/)(31-2/) 


• ( 68 ) 


The same equations can also be used to obtain equa¬ 
tions for the RG invariant ah (also announced in |12|). 
Eqs. ED, ED and 1631) imply the NLO formula 


NLO 


2 - (33 -2/) In ^ 

07r [i 


153-19/ 

2tt(33-2/) 


In 


l + ^(33-2/)ln^ 

07T fl 


(69) 


Again, Eq. ED can be used to write ln(m/j/ p) in terms of 
\n(mh/p). This leads to the following asymptotic formula 
for ah as mh —» oo: 

ah ~ j |(33 — 2/) In + Kf lnln + 0(1) j, 

3(153- 19/) 12(33- 2/) 

Kf = 33-2/-31-2/ ■ (70) 

These results show that we have complete control over 
the matching process. Once closed expressions such as 
ED and ED have been obtained, RG invariance can be 
maintained for each logarithmic order, and so there is no 
need to truncate to a given order of perturbation theory. 


VII. MASS-MATCHING FUNCTION 

Most of the analysis above is restricted to the case of 
just one heavy quark h , but it can be readily extended 
to include sequential decoupling, where heavy quarks are 
decoupled one at a time. The new feature which arises is 
the need to match the mass of the second heavy quark. 
For example, suppose that, having decoupled the t quark 
in F = 6 flavor QCD, we would like to decouple the b 
quark as well: 


rrit —> oo first, then mb —> oo. (71) 

Then it will be necessary to match the six-flavor defi¬ 
nition Wbe = mb of the bottom quark mass to its five- 
flavour definition mb§. 

As for the matching of couplings, the key is to find a 
RG invariant definition of mass to which the Appelquist- 
Carrazone theorem { 6 ] can be applied. This problem was 
solved by Bernreuther 0, again by recourse to the mo¬ 
mentum subtraction (MO) scheme. 

Let Aiip 2 ) and Be(p 2 ) denote the form factors for 
the unrenormalized 1PI light-quark self-energy amplitude 


—i(jfiAe — moeBc). This corresponds to the unrenormal¬ 
ized quark propagator 


S ^ $1(1 -A t ) - m oe (l-Bi) 

= ^-mo t (l-Bi)/(l-Aty (72) 

Define MO light-quark masses Mf AO (Q) at a fixed space¬ 
like point p 2 = — Q 2 : 

My io {Q) = m oe (l-Bi(-Q 2 )) / (l-Ai(-Q 2 )). (73) 

This mass is RG invariant because, expressed in terms of 
renormalized quantities, it is finite. The choice of space¬ 
like subtraction point — Q 2 means that mass renormal¬ 
ization, as well as coupling-constant and wave-function 
renormalization, is performed off-shell : 2 

iS~ 1 (p )\ p2= _ Q2 =i>~ Mf°{Q) ± 0 at p 2 = -Q 2 . (74) 

This avoids problems with Bloch-Nordsieck logarithms 
produced by n-loop perturbation theory at the on-shell 
point p 2 ~ nr 6 |: 


iS 1 (p)~(p 2 —m 2 ) 1 [ln ($» 2 — m 2 )]". 

A complication familiar to many authors (T71 ) is that, 
unlike an on-shell renormalized mass, M^ mo (( 3 ) is gauge 
dependent. Despite this, the resulting mass-matching re¬ 
lation between mef and mpf+\) is gauge invariant 3 
because the relations between MS masses and their bare 
counterparts are gauge invariant. In two-loop perturba¬ 
tion theory, the result 0 for QCD with three colors and 
/ light flavours is 


me/ 
m i(f+ 1 ) 


= 1 


a 


/+! 

12tt 2 


In 2 


p 2 


^ln^ 

3 p 2 


89 

36 


+ °( a /+i)- 


(75) 


We would like to extend this result to include all terms 
of the same logarithmic order. This is achieved by intro¬ 
ducing our second matching function, Q - the matching 
function for mass renormalization. On order to reduce 
notational complexity, we consider the special case / = 5 
mentioned at the beginning of this section. 

Consider the RG equation 


V 6 ln^ = 6 6 (ae) - S 5 (a 5 ) (76) 

m b5 

The leading power in a large -mt expansion of mb 6 /mb 5 
is a function of 0:5 and In (mt/p) but does not depend 


2 


A similar MO definition for heavy-quark masses M ^ E3 yields 
a /3-function (3(g , Mh/Q) El with smooth threshold behavior at 


Q ~ M h . 


















on light-quark masses, so the general solution of Hi is 

in^ = r^M-f^TTS 

m b5 Ja 5 P5[X) J as p 6 (x) 

+ g 6 ^ 5 (a t ) + 0(rn t _1 ) (77) 

where t/6->s(5t) is the mass-matching function. Like 
iF 6 ^ 5 ( 5 t ) in Eq. 14511 . t/ 6^5 arises as an integration con¬ 
stant of a RG equation, so it can depend only on the RG 
invariant a t . Also like T , it cannot be deduced from the 
RG and must be calculated separately. 

At one-loop order, there are no top-quark corrections 
(Fig. ID, so the mass-matching condition is trivial 0 : 

m mo = + a 6 {l-loop}] +0(al) 

= m^s [1 + a 5 { 1-loop}] + 0(a 2 ). (78) 

Here a/{1-loop} denotes the self-energy amplitude de¬ 
rived from Fig. □ Eliminating M b MO , we find 

l n mb 6 _ ^ _ a 6 ){i_i 00p } + 0(al) = 0(a§). (79) 
m b5 

At this point, we need to specify what is LO, NLO, 
and so on. If we were talking only about corrections 
to mass, we might consider terms ~ (aflnmt/p) as 
LO, but in general, it is more convenient to regard them 
as NLO. That is because mass renormalization does not 
contribute to physical amplitudes in LO. With this con¬ 
vention, Eqs. CHI . iTvSli . and iTTfll imply 



s f(x) 

Pf(x) 


NLO 


NLO 


6 

33 

2 J 


6 


33 

2 






M . 
(80) 


In lowest-order perturbation theory, we have 

Ja 5 Psix) J ae f3e{x) 

= ^(a 6 -a 5 )ln^+0(ag) = 0(af) (81) 

since the couplings and 0:5 differ by 0(a|). From 
Eqs. (1771 . (1711 and EQ), we conclude 

06^5 (on) = 0(a%). (82) 


Having established that Q is irrelevant at NLO, we 
neglect it in Eq. o and substitute This yields the 


b 


complete NLO expression: 


In 


m b 6 _ 6 

m b5 nTo 


In 


1 a6[33_ 

+ 3tt ^ 2 


6 ) In 


In 1 4- —— 


a b f 33 


m t 


m t 


If desired, ag can be eliminated in favour of a .5 via 
Eq. (13511 . Note that the 0(a|) NLO term is a double 
logarithm which arises from the diagram with a f-loop 
inserted in the gluon propagator (Fig. [21 : 


ln m b6 _ 2 

77165 NLO 7r 

+ 


(ag — a,) ln — 3^2 (" 2 " 

^5 ( 33 -A , 2 rrit , 3 , 

sj 111 y + °(“=) 


NLO 


I GW 

3 V 7T / fj. 


0(al). 


ln 2 ^ 


(84) 


This reproduces the NLO term of Bernreuther’s result 
Gil for 3 rrif/mf+i. 

Eq. m generates the complete set of NLO logarithms. 
They correspond to diagrams with a string of one-loop 
f-quark bubbles inserted in the gluon propagator of the 
one-loop 6 -quark self-energy amplitude. 

The constant term 89a|/(4327r 2 ) in Eq. itTbl) corre¬ 
sponds to the first non-zero contribution to the matching 
function Q. We state the result for any value of /: 

oq 

<7 /+w (5 h ) = + 0(5 3 ). (85) 

This term is required if NNNLO corrections are being 
calculated. 


VIII. APPLICATION TO SEQUENTIAL 
DECOUPLING 

When decoupling the b quark, it is natural to define 
five-flavor quantities a b5 and fh b5 by analogy with the 
six-flavor running coupling a t and mass frit f° r the top 
quark: 


t 



b 


FIG. 2: Two-loop lieavy-quark contribution to mass renor¬ 
malization of the bottom quark. 


3 Note that Bernreuther expands in ae instead of as, but to this 
order the coefficient is the same. 


FIG. 1: One-loop correction to mass renormalization. 
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In mb5 

r&b 5 

/ dx 

J CK 5 

1 - d 5 (x) 

( 86 ) 

ft 

feix) 

In 77165 - 

/•“t 5 

/ dx 

h(x) 

(87) 

m b5 

J 015 

p5{x) 


Clearly, both a b5 and m b5 are RG/ = s invariant, 

T> 5 a b5 = 0 , T> 5 fh b5 = 0 (88) 

but in fact, they are also RG / = 6 invariant as a conse¬ 
quence of Eqs. O and d: ' 


V 6 a b5 = 0 , V 6 m b 5 = 0. (89) 

This means that a b5 and m b5 can be expressed in terms 
of invariants of the original six-flavor theory . 4 To see this, 
first combine Eqs. El and <2D> and then CD: 


m t 


m b 5 


In —- — In —— = / dx 

mt m b5 


6 6 (x) 


(k(x) 

= ^ln^ + 

m b5 J otb 5 




/%(a 


ix MT) 


+ Ge^5(at) + 0(m t x ). 
This equation simplifies to 


m t 


In = In ( 

m b 5 


— ) + 

\m b J 6 . 


at S 6 (x) 

Ps(x) 


OLb 5 


+ f?6^5(cq) + 0 (m t *)• 
The sum of Eqs. and 1871) 


In 


m b5 _ f 

P J a 


a b5 l 

a 5 dX k{x) 


can be subtracted from Eq. ©, with the result 


dx 


m b5 


'<*b5 


Pb{x) 


(90) 


(91) 


(92) 


+ f 6 ^ 5 (a t ) + 0{mt l ). (93) 


If (1931) is now combined with ED- we find that a b5 and 
hence fh b5 can be expressed in terms of RG invariants 
of the original six-flavor theory, viz. at and the ratio 

(■ m t /m b ) 6 : 


In(^) = 

\m b J 6 Ja b5 As 0*0 

+ ^ 6 ^ 5 ( 5 () - £ 6 —5(a t ) + 0(m t _1 ). (94) 


4 This property is essential for any generalization of the analy¬ 
sis to simultaneous decoupling, where all couplings and masses, 
running or otherwise, will have to be defined only in terms of 
the initial theory (F = 6) or the residual theory (/ = 4 if just 
the t and b are being decoupled), with no reference to five-flavor 
couplings or masses. 


The sequential decoupling of the t and b quarks refers 
to the limiting procedure 

ln(m t //i) 3> \n(m b /p) —> oo (95) 

where we choose a six-flavor definition for m b as well as 
rri t . Leading-power six-flavor amplitudes are represented 
by logarithmic expansions for t-quark decoupling 

•Ae ~ y~] C p5 ln p (m t /p) (96) 

p>0 


where each five-flavor coefficient C p 5 is a leading-power 
asymptotic expansion for 6 -quark decoupling: 

Cp 5 ~ 1 n q (rn b /ji) (97) 

q^Q 


The last decoupling (that of the b quark) is carried out 
with a 4 held fixed. Therefore we seek formulas for cou¬ 
plings such as «6 and at in terms of mt, m b and 0:4. 

As always, the key step in the derivation of NLO for¬ 
mulas is the neglect of some matching functions. In this 
case, we neglect the NNLO functions T, 6^5 = 0(a t ) and 
4 = 0(a b5 ) for matching a 6 to a§ and a b to a 4 , 
and the NNNLO function ^ 6 _>5 = 0(af) for m b5 to be 
matched to m be = m b . 

We start with the NLO formula (El for the five-flavor 
mass m b5 . To this order, all dependence on <25 and ae 
can be eliminated via LO formulas derived from Eq. El, 


a 5 = a 4 


/( 


1 + — In —) 
37t h J 


(98) 


“ s it “ 4 /( 1 + S : ( ln y +1 "y)} <"> 

where (again to this order) m b5 may be replaced by m b 
on the right-hand side. Then Eq. El becomes 


In 


m b5 


m b 


12 


H NLO 


= In - = - — — In (1 + 


a 4 m b ‘ 




23 


37T 


■In 




161 \ 67t fj, 3 tt J 


+ -j In ( 1 + 


o ?4 m b 


37T 


■In 


a 4 


m t \ 




+ 7T“ In —] • (100) 

07T fl J 


Similarly, consider the NLO relation El between 
a /+1 and a/. For the case f = 4, the heavy-quark mass 
TOh in (1681) is m b5 , but we can use Eq. (HOOD to eliminate 
m b 5 in favor of m b , 


a 5 1 = a, 1 
b NLO 


1 ln TO fc5 
37T p 


+ C4 In 


NLO 


a 4 m b 

l + u ln ~ 


d 4 In 


25a 4 . m b 

1 + —— In — 
07r fi J 


/i J 

( 101 ) 


where the constants c 4 and d 4 are given by 


c 4 = 33/(23tt) , d 4 = 121/(1150tt). (102) 



















For / = 5, nth in (PI) is the six-flavor mass to*. Any 
dependence can be removed via Eqs. or iron : 
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Ckc — CXc, 
b NLO 


-1 


NLO 


1 i TOt , i 

— in — + C 5 In 

07T H 




LO 


+ d 5 In 


1 + 


23a 5 


67 r 


LO 


, rn t 

In — 
/' . 


, rn t 

in —— 

(103) 


The coefficients C 5 and d b have numerical values 
c 5 = 47/(42tt) , d 5 = 137/(966tt). 

The result of combining these formulas is: 


(104) 


a 6 1 = «4 

0 NLO 


1 + — (in — +ln — 
37 t V ft n 




55 
427T 
121 
9667T 
121 
11507T 


1 / 1 1 1 Tfib 1 wit 

m ( I + 7 — In — H-In — 

37T fj, 37T fl 

, 23a 4 , m t a 4 m b \ 

in I 1 + —;— In — + — In — 
07T fl T fi J 


25a 4 . rn b 

ln ( 1 + — 7 —In — 

D7T /i 


(105) 


The same procedure can be applied to the NLO for¬ 
mula m for the RG invariant mass nih- For / = 4, the 
result is 


In 


m b 5 


/j, NLO 


= ln^-Sta(l + ^lln» 


m b 


12 

23 


67T 


161 


In [ 1 + In — 

07T /i 


OL\ 


m b \ 
3n £l ) 


In 


+ -j in ( 1 + 


aq m b 


37r 


■In 




a 4 m t 

— in 
) 1 /: n 


(106) 


after substitution of Eq. cm for m b 5 . For / = 5, Eq. 
<E3 expresses mt in terms of 0 : 5 , for which the LO for¬ 
mula (tUSl) may be substituted, with the result 


m t . m t 4 / a 4 . m t a 4 . m b 

in = in —— -|— In 1 H-In ——|-in — 

l_i NLO i_i 7 V 37T /i 37T fi 


,»( 


1 + 


/i 07T /i 

23a 4 mt a 4 m b \ 

-In ——|-In 

1 :■ ' fi 37r fi / 


• (107) 


If desired, inverses of m and m can be constructed 
and used to express quantities such as a§ in terms of the 
invariant masses fh bb and fh t instead of m b and m*. 

Finally, we extract NLO formulas for the invariant run¬ 
ning couplings for sequential t, b decoupling from Eq. PI 
for ah- For the case f = 4, Eqs. lIHTfl) and (HOOD imply 


-1 


cthr. ± = a, 1 H-In —- — 

5 NLO 4 07 T fi 115077 


25 m b 


729 


In 1 + 


25ct4 m b 

-In 

1 5" p, 


100 

483tt 


In 1 + 


23a 4 TO* a 4 m b 

— r — In — + — In- 

1)77 fi 3n fi 


+ ^-ln( 1 + |lln^ + filn^). (108) 

217T \ o7T fl o7T fi 1 


For / = 5, it is necessary to combine Eq. © with Eqs. 

front , (trnri) and (1TU7I : 


= a 4 

NLO 


23 m t 1 , m b 
— In — + — In — 
07T /i 37T fi 


457 



483tt " 


50 , 


21 tt 

+ 

121 


In 1 + 


23a 4 mt 

— r —In — 

07T [l 


(%4 i Wl b 

X- In — 
37T fi 


1150tt 


ry 4 , m b a 4 mt 

1 -|-In —-|-in 


, 25a 4 , m b 

in ( 1 + ——In — 

07T zi 


(109) 


A useful check of the formalism can be obtained by show¬ 
ing that the difference 


a t a b5 

is correctly given by Eq. (F7H) in NLO. 


IX. THEORIES LACKING ASYMPTOTIC 
FREEDOM 

So far, we have limited the discussion to heavy fermions 
in a gauge theory such as QCD and used asymptotic free¬ 
dom to obtain decoupling in the infinite-mass limit. 

Does decoupling work if asymptotic freedom is not 
present, as in quantum electrodynamics (QED) or scalar 
field theory with X(f > 4 interaction in four dimensions? Per¬ 
haps this is part of the wider debate [^3 about whether 
such theories are inconsistent or trivial, particularly for 
the continuum limit of the lattice approximation. 

We start from the premise that a theory makes sense 
only if its heavy particles decouple in the infinite-mass 
limit to produce another consistent theory. Questions 
about how the non-perturbative theory could depend on 
details of regulators and their removal are not consid¬ 
ered. We simply assume that some means of produc¬ 
ing a fully interacting cutoff-independent theory has been 
found, e.g. for QED& with k species of equal-charge lep¬ 
tons, and apply our premise. 

The notation is similar to that used above for QCD. 
Let CKfc = e|/( 47 r) be the MS renormalized fine structure 
constant for QED^, where the charged leptons have MS 
masses mj, j = 1 ,..., k, and let ( 3 k and 5 k be the Callan- 
Symanzik functions for charge and mass renormalization. 
Denote by an and fha Witten’s invariant versions of the 
running fine-structure constant and heavy-lepton mass. 
Then, repeating the arguments of Secs. rm and Em we 
can construct coupling-constant and mass matching func¬ 
tions T and Q for the decoupling of one species of lepton: 


QED fe+1 —> QED fc , k ^ 1 (HO) 

The free-photon theory QEDo lacks a /3-function so it is 
a special case. 
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First consider the analogue of Eq. m, including the 
case k = 0: 


X. CONCLUDING REMARKS 


% 

P 


r-a H 

/ dx 

' aic+l 


1 

Pk+i(x) 


The decoupling condition 


( 111 ) 


otk+i —> 0 as mu —> oo for fixed (112) 

involves the x = 0 solution of the equation j3k+i{x) = 0, 
but it produces an infinity of the wrong sign because this 
fixed point is infra-red stable. For consistency, olh must 
approach a singularity of the integral sufficient to reverse 
the effect of the x = 0 contribution. This could arise 
from an ultra-violet fixed point, 


&H * &k-\-l,oo (113) 

or else q.h approaches x = oo, in which case we must 
suppose that l//3k+i{x) is not integrable. 

Next define T as in Eq. ©, with k ^ 1: 


The introduction of the matching functions T and Q 
for coupling constants and masses (Eqs. © and CZ3) 
completes the theoretical structure needed for a system¬ 
atic application of the RG to the decoupling of heavy 
particles. We have considered just QCD and QED, but 
the field-theoretic principles are much the same for any 
theory. The main case still to be checked is a full RG 
analysis of the decoupling of heavy gauge bosons whose 
masses are induced by the Higgs mechanism. 

In this article, we decoupled only one particle at a 
time (for simplicity) and concentrated on field-theoretic 
aspects of the subject. Actually, our work on matching 
functions arises from a need to consider the simultaneous 
decoupling of more than one heavy quark in phenomeno¬ 
logical applications. The challenge is to keep track of 
dependence on large logarithms such as 


, rn t . m t , m b 

In- = In —-in — 

m b \i /j, 


(116) 


fflH 1 

ln^—= / dx—— + T k+ i^ k {a H ) (114) 
P J ak Pk(x) 

In the decoupling limit irm the singularity on the left- 
hand side is generated entirely from the running of an- 
We can conclude that this is due to a QED fc fixed point 
ctk ,oo only if it coincides with that of QEDfc + i: 

l,oo — CKfc, oo (fffi) 

Similarly, a non-integrable singularity at x = oo can be 
the sole cause only if this happens for both QED*, and 
QEDfc + i. Otherwise, we must suppose that the matching 
function has a singularity at afe j00 or oo. 

The non-perturbative theory of QED of Johnson, 
Baker, Willey and Adler 0j is an example of the case 
(111511 . Indeed, if the arguments of Adler for an infinite- 
order zero at the fixed point are applied to a many-species 
theory (all with the same charge), the result is clear: 
there is no dependence on the number of species. 

Notice that these conclusions are driven by the lack 
of asymptotic freedom of the initial theory. For exam¬ 
ple, consider what happens in QCD when a heavy quark 
decouples from the non-asymptotically free 17-flavor the¬ 
ory to produce the 16-flavor theory with asymptotic free¬ 
dom. Eq. © for F = 17 implies that a b increases. Thus 
in Eq. © with / = 16, ah is driven towards a non- 
perturbative infra-red region of the residual theory. 


Since these logarithms do not depend on the MS scale p, 
the conventional tactic used in phenomenology for sin¬ 
gle heavy quarks fails: there is no way of making such 
logarithms small by choosing p = 0{nrih)- 

As indicated in our work jl 2j on NLO heavy-quark ef¬ 
fects in axial charges of nucleons, the analysis can be 
generalized to include simultaneous decoupling of several 
heavy particles. This includes the introduction of match¬ 
ing functions of several variables, one for each heavy par¬ 
ticle. We will present this extension of the theory in a 
forthcoming publication. 

One can also anticipate generalizations to situations 
where momentum and mass logarithms grow large to¬ 
gether. Examples from the literature occur in collider 
physics , Higgs and supersymmetric particle produc¬ 
tion 0|, and deep-inelastic scattering through thresh¬ 
olds for heavy-particle production 0}. 
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